Present status of the theory of the high T c cuprates 
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The Gutzwiller projected mean field theory, also called Plain Vanilla or RMFT, is explained and 
its successes and possible extensions in describing the phenomenology of the cuprate superconductors 
are discussed. Throughout, we emphasize that while this is a Hartree Fock based BCS theory, it 
embodies fundamental differences from conventional perturbative many body theory which may be 
characterized by calling it a theory of the doped Mott insulator. 



I. HISTORICAL NOTE 

In early 1987, just as the remarkable Bednorz-Muller 
discovery was becoming widely known, the basis for the 
theory of the materials which they had discovered was 
laid downi. It was observed that the Cu02 planes on 
which they are based were plausibly describable by a 
particularly simple version of the Hubbard model, the 
case of a single non-degenerate band, and that the "sto- 
ichiometric" case where the nominal valence is Cu ++ is 
well described as a Mott insulator. The superconduc- 
tors are obtained when, in the "reservoir" layers between 
the planes, substitutional impurities of lower valence are 
introduced, thus doping extra holes into the Cu d-shell 
(which is of course strongly hybridized with the O p- 
shells, according to the well-known principles of ligand 
field theory.) All of the plausible theories about these 
materials describe them as "doped Mott insulators". 

A mechanism for electron pairing in mixed valence 
systems, which are somewhat similar, already had been 
suggested by two groups^, namely using the antiferro- 
magnetic "superexchange" interaction between spins as 
a pairing force. In Ref. 1, I likened this pairing force to 
the valence bonding effect for which it is essentially re- 
sponsible, and pointed out that the old idea of a quantum 
liquid of valence bonds resonating around among differ- 
ent pairings of atoms had a great similarity to supercon- 
ductivity. In fact, I proposed an explicit form for such a 
state in terms of a Gutzwiller-projected BCS paired wave 
function, and in a series of papers in 1987 elaborated on 
formalisms for getting continuously from the Mott insu- 
lator to the superconductor. 

Unfortunately, through a series of misjudgments on my 
part, which are permanently recorded in an unfortunately 
timed bookS, my group and I thereupon fell off the cor- 
rect trail to a solution, only to return to the correct path 
ten years later once we had absorbed the unequivocal 
experimental evidence that my "interlayer" theory was 
wrong. But fortunately, at least two separate groups had 
in the meantime built a theory on the 1987 foundations 
which turned out to be basically corrected. In this article 
I will follow the second of these references but they are 
equivalent. The important thing about both is that they 
realized that the correct solution of the original undoped 
RVB problem was not the isotropic "extended s" which 
I had been discussing but a more complex one with both 



s-like and d-like gaps, which Kotliar called u s + id" . Both 
of these papers predicted the real d-wave gap with nodes 
which was eventually observed, and in addition a number 
of other results which were to be confirmed one by one in 
the coming years. It has been our perverse fate that the 
theory, properly handled, has made one after another cor- 
rect prediction, well ahead of the experiments, but that 
these have been obscured by irrelevancies and misinter- 
pretations until the mistaken impression has arisen that 
the whole subject is utterly mysterious. 

It was not for another 5 years that the d-wave gap was 
verified, and by that time the field had suffered from a 
proliferation of proposed theories of greater or lesser de- 
grees of plausibility. The gradual experimental unveiling 
of the facts about the cuprates sometimes meant that 
each experiment came with a built-in theory and that 
theories which had predicted the result long before were 
not sufficiently "up-to-date" to enter the public discus- 
sion. For instance, the d-wave came to be identified with 
the idea of propagating "antifcrromagnetic spin fluctua- 
tions", which was a popular fad at the time of its ver- 
ification, rather than with its earliest, and much more 
natural, prediction in Ref. 00. Another example of this 
phenomenon was the observation of the "spin gap" or 
"pseudogap" in underdoped materials above the super- 
conducting " dome" , again an obvious consequence in Ref. 

but as it revealed itself it received a congeries of 
faddish explanations from local theorists: a mysterious 
"quantum critical point" , a "spin nematic" , again AF 
spin fluctuations, the "d-density wave" , you name it. 

In any case, these early theories only came to be re- 
vived in the early 2000 's by groups which were able to use 
them as the basis for accurate quantum Monte Carlo cal- 
culations using realistic parameter values^, and brought 
forward without too much modification some of the pre- 
dictions which had looked so surprising in 1988 but had 
been very close to correct. A group of us summarized 
the successes of the theory, adding a small amount of 
further physical ideas, in a review paper which we called 
the "Plain Vanilla" theory of high T^. Here I will review 
that theory and the subsequent developments, including 
particularly the explanation and calculation of asymmet- 
ric tunneling spectra using it, and the recent theory of 
the pseudogap phase which throws a great deal of light 
on the overall physics of the phase diagram. 

Since a great deal of emphasis has been put on the 
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problem of the epistemics of complex phases like the high 
Tc cuprates, and whether a meaningful solution to the 
accompanying puzzles can be found , I'd like to spend 
a few sentences on that aspect. First, a bit about the 
nature of condensed matter physics. Among the sciences 
this one is almost uniquely ovcrdctcrmined, experimen- 
tally because of the variety and precision of the probes 
which can be applied, and theoretically because the quan- 
tum physics of atoms and electrons is so well understood. 
I have always maintained that the correctness of a theo- 
retical hypothesis is assured in this field if it can find a 
way to fit in with all these constraints: that there is likely 
to be only one possible way to fit all-or even a majority-of 
the observations together, and not to violate any theo- 
retical impossibilities. In this process of fitting things 
together there is no room for one-experiment theories, 
doctrinal conservatism (the older generation and some 
younger scientists won't let go of phonons), or yet un- 
trammeled imagination (anyon superconductivity, SO (5), 
QCP's, perhaps interlayer tunneling). The naked reality 
is strange enough. 

A final word. The way you know you are right is when 
you wake up and realize that you have the answers to 
deep, fundamental questions that you didn't really know 
to ask or expect to answer. For the old superconduc- 
tors, such a question was "why are polyelectronic metals 
favored?" -the question Pines, Morel, myself and McMil- 
lan answered with dynamic screening for the phonon the- 
ory. Here there are at least two such questions: "Why 
the cuprates-what is unique about copper?" ; and, "Why 
d-wave and why is the gap persistently real? That is, 
why the striking nodes?" The second is the question I 
didn't think to ask, but it is profound-any other simple 
mechanism which leads to a d-wave can lead also to an 
xy or isotropic symmetry, which will appear in quadra- 
ture in order to fill in the nodes, which are intrinsically 
unstable in a BCS theory. The mechanism by which the 
A phase of 3 He acquires nodes was, for instance, crucial 
to our understanding of that system. 

II. THE PLAIN VANILLA (RMFT) THEORY 

The underlying concept of the plain vanilla theory 
is very simple. In fact, it follows as closely as possi- 
ble the precedent of the BCS theory. The BCS the- 
ory in its original form is a generalization of Hartree- 
Fock theory to allow for not only the direct and ex- 
change mean fields, which appear in the one-electron 
mean field Hamiltonian as V(r)p(rO = v^*(r)^>(r), and 
A(r, r')^f*(r')^f(r), but also the "anomalous" self-energy, 
A(r, r')\E , *(r)5'*(r') + h.c.. These result from the three 
possible ways to factorize the interaction energy, 

J V(r-r')^*(r)^(r)^*(r')^(r')drdr' . 

BCS theory is basically a variational theory: the assumed 
wave function is a simple product of one-quasiparticle 



operators creating quasiparticles from the vacuum, and 
the "gap" equations, equivalently to the mean field equa- 
tions, determine that the quasiparticle creation opera- 
tors all have positive energies, so that all possible single- 
particle excitations increase the energy. 

Simple Hartree theory won't work for a Hubbard model 
in which the on-site interaction energy is the largest en- 
ergy in the problem. Very early on£, it was realized that 
the solution to that problem was to transform to a rep- 
resentation in which the on-site interaction energy U has 
been renormalized to oo as opposed to the conventional 
scheme well described by Shankar 9 where the idea is to 
transform to some system of noninteracting entities. We 
employ a canonical transform expiS to eliminate all ma- 
trix elements of the Hamiltonian which lead into the sub- 
space in which two electrons simultaneously occupy the 
same site; i.e., those which have the large energy U. This 
transformation 8 can be derived perturbatively as a series 
in inverse powers of U. That is, we start from the "real" 
Hamiltonian 

H = ^2 ^i c L c jo- + n - c - + u ^2 n i\ n il > (!) 

along with direct exchange and smaller terms, and trans- 
form it into the t — J Hamiltonian; Hq — ► e lS Hoe~ lS — 
H t -j + 0(t 3 /U 2 ) + . .., where H t -j is given by, 

H t -j = ^ tjjPGc\<r C j<r P G+h.C.+y^ j Jjj (Sj-Sj- jUjllj) . 

(ij) 

In the above equation, we have ignored terms includ- 
ing longer-range Coulomb and phonon interactions-which 
latter are not particularly small, but clearly are incapable 
of causing the gigantic superconducting gaps which are 
observed. Here, Pq is the full Gutzwiller projector which 
hereafter we will call P: 

p = na - • ( 3 ) 

i 

That J2J) is really a correct description of the electronics 
of the cuprates was tested first by Schluter et alJ&, in 
1988, who found that the calculated energies of low-lying 
states in small clusters of the cuprate structure, using 
the full Hamiltonian, were well reproduced by the trun- 
cations implicit in equations (JIJ and (J2J). (Another early 
discovery long since forgotten.) 

Always remembering that the t — J model wave func- 
tion must be transformed by exp iS in the end to repre- 
sent Hubbard model reality, we proceed to try to find a 
variational ground state for Clearly, since the Hamil- 
tonian is now in block diagonal form, any low-energy 
state must contain only amplitudes for the projected sub- 
space, so that, 

V = P$(r 1 ,r 2 ,...,r N ) , (4) 

where $ is a general TV-particle wave function. The 
essence of the "plain vanilla" approximation is to pro- 
pose that we approximate the wave function to be 
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projected, using the Hartree-Fock-BCS ansatz that it is 
a product of quasiparticles. I can see no reason that this 
is apriori less reasonable than the BCS theory itself. If 
there is a single-particle-like representation of the ground 
state, this is the way to derive one. In the event, there is 
such a representation, experimentally-by now there is all 
kinds of evidence that the state has gapped quasiparti- 
cles near a large Fermi surface, over a fairly wide range of 
doping-the so-called "dome" region of the phase diagram 
of T vs doping. I can't too much emphasize this: this 
procedure is the natural, and probably the only, way to 
derive a BCS-like superconductor from the t — J Hamil- 
tonian. 

A second, and less certain, fact is that the resulting 
excitations may be reasonably sharp and well-dcfincd- 
though, because of the projection operator, the same 
may not be said of actual quasiparticles : cP is not the 
same as the single-particle-like excitation Pc. But the 
representation in terms of Pc's has some as yet unre- 
solved peculiarities: it is overcomplete, which may mean, 
among other things, that the excitations can scatter each 
other very strongly. But the fact of overcompleteness 
does not much affect either the variational equations nor 
the validity of them as giving the energies of approximate 
single-particle excitations. In writing out these equations 
we follow Ref. in self-consistently choosing a particu- 
lar relative gauged for the J-term relative to the kinetic 
energy. This choice is discussed later. 



the Renormalized Mean Field Theory: 



Our Ansatz for <I> in Eq. J2J ,then, is 
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£,k = gtk + <a = .9£fe + gjJ ^ 7fe - fc ' 
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(6) 



Here, Jjk-k' is the Fourier transform of the exchange 
interaction (assumed nearest neighbor) ek is the bare, 
unrenormalized kinetic energy, A and £ are the anoma- 
lous and normal self energies, £, the renormalized kinetic 
energy and Ek is the quasiparticle energy. 

In Fig. ^ we present results for the magnitude of the 
d-wave gap, A, and the size of the order parameter from 
Edegger et al.'s solutions of the gap equations 12 , just to 
convince the reader that these track the observed maxi- 
mum gap and dome reasonably well. A generalized phase 
diagram incorporating the results of a number of experi- 
ments is shown in Fig. [21 (This figure differs from a phase 
diagram often drawn for which the T* line intersects the 
dome and no trace of the pseudogap phase remains for op- 
timally doped materials. Ong's Nernst effect data among 
others seem to unequivocally reject this interpretation.) 
Since 1988, it seems, the quantitative explanation of high 
T c superconductivity has been available. 



III. EXTENSIONS OF THE RMFT 



(5) 



A. Spin-Charge Locking 



In the recent papers by Paramekanti et ai, the param- 
eters u and v were evaluated variationally using Vari- 
ational Monte Carlo techniques^. But the results were 
almost identical to those found in the earlier papers us- 
ing a very simple approximation due to Gutzwiller, which 
is exact in the limit that the gap is small relative to the 
Fermi energy. In this approximation we assume that the 
correction to the probability of occupation of the sites 
caused by projection is uncorrelated spatially, because, 
obviously, the projection operates only site by site, ig- 
noring the occupancy of neighbors. Thus the correction 
may be estimated by simply calculating what happens 
to the average occupancies. It is easily shown that the 
change in the average number of neighbors with one site 
empty, the other singly-occupied, is a reduction by the 
factor g = 2a;/ (1 + x), while the change in the num- 
ber of pairs of singly-occupied sites is an increase by 
gj = 4/(1 + xf = (2- g) 2 . Thus the effect of the kinetic 
energy is reduced by the factor g, and that of J is in- 
creased by g J , but otherwise, in this approximation, we 
employ the t — J Hamiltonian in precisely the same 
way as a real one. Thus we arrive at the "plain vanilla" 
gap equations in the "Gutzwiller approximation", i.e., 



Note that as g — > (the "true" RVB), A and <; are 
interchangeable. This represents a deep reality: that 
for the half-filled Mott insulator, the representation of 
the magnetic state of the spins by fermionic variables - 
the "spinons" of RVB theory-is doubly overcomplete . 
One may represent an | spin on site i either by creat- 
ing an | spin, cj^ , or by destroying a j on that site, 
Cjj_, or by any unitary superposition of the two. In terms 
of a hypothecated RVB state, described as a Gutzwiller 
projected BCS wave function at half filling, this means 




FIG. 1: (a) Doping dependence of the dimensionless mean 
field parameters f, A; (b) Doping dependence of (solid) 
the SC order parameter, <3>, and (dashed) the gap, |Aj.|, at 
k — (-7T, 0) in units of t. 
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that the three Anderson-Nambu spinors r, (i — 1,2,3) 
of the BCS state may be rotated at will, since they rep- 
resent quantities which transform into each other when 
the SU(2) transformation is applied. The constraint of 
the Gutzwiller projection also requires that only two of 
the three r vectors have finite self-energies attached to 
them, so that the symmetry is fully expressed as local ro- 
tation of a dyad of self-energies A and £ which must be 
perpendicular to each other. All of the various alterna- 
tive states which have been proposed-the "flux phase", 
the d-density wave, the staggered flux phase, etc., are 
one or another of these totally equivalent states, in the 
half-filled case. The two "gaps" , for the minimum-energy 
solution of Eq. |JB}, are of maximally different symme- 
tries. In the half-filled case, and in the special case that 
we have only nearest neighbor exchange so that 7 is of 
the form cos k x + cos k y , the two are equal in magnitude 
and of the form cosk x ± cosk y . The only point where 
both vanish is where both fc's are 7r/2, which gives the 
nodes which are the common feature of all the equivalent 
"ghost states" I mentioned above. 

It is irrelevant that the actual half-filled band is not the 
RVB state but a commensurate antiferromagnet, which 
has slightly lower energy for the Heisenberg model. It 
is still meaningful to examine the solutions of the full 
gap equation by referring them back to the hypothetical 
limit g — 0. What happens is that, as we reintroduce 
the kinetic energy by doping, the antiferromagnetic state 
does not gain kinetic energy as rapidly as the best RVB 
state, and the latter prevails at a few percent doping. 
Actually, the equations © represent a special choice of 
gauge, and we could in principle orient the kinetic energy 
along any chosen axis in the r-space, and minimize the 
energy as a functional of that orientation-the resulting 
equations are given elsewhereii. But it is clear that the 
optimum kinetic energy is achieved when the "£" axis, 
the function with the symmetry cos k x + cos k y , is chosen 
as an ordinary self-energy as in © . Then the other form 
of solution, the odd combination cos k x — cos ky , act s m 
the direction t\ and serves as an anomalous self-energy or 
gap function. This is the principle I called "charge-spin 
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FIG. 2: Generalized phase diagram of the cuprate supercon- 
ductors 



locking" 11 . The locking energy was estimated in that 
reference as well as by Kotliar and Liu and found to be 
large: of order gt for small dopings and comparable with 
T* for larger ones. 

This large locking energy means that the gap structure 
is established at temperatures well above the supercon- 
ducting "dome" of T c 's. The reason the system does 
not become superconducting is that the phase stiffness 
is weaker, at least for doping up to the optimum, than 
the gap energy, in contrast to the BCS case. T c is deter- 
mined by the proliferation of vortices, not by the break- 
down of pairing. Experimentally, in systems which are 
basically two-dimensional, one sees Kosterlitz-Thouless 
transitions; and the cleanest measurements for optimal 
YBCO find 3D XY model exponents, very accurately^. 
Both observations indicate that the order parameter am- 
plitude remains finite above T c , and in fact the observa- 
tions of Ong on Nernst effect and nonlinear diamagnetic 
susceptibility^ show that a vortex liquid state persists 
well above the dome, especially on the underdoped side. 
From these measurements, as well as theory, we are be- 
ginning to establish that what has been called the " vortex 
liquid", i.e., a disordered superconductor as opposed to 
a normal metal, may be a distinct state of matter which 
is particularly characteristic of the cuprates. 

That T c embodies a transition to a vortex liquid state 
suggests a phenomenology of this metallic state above 
T c quite different from that of a normal metal. We must 
think of it as everywhere superconducting, but filled with 
a tangle of thermally-generated vortices (at low fields.) 
The supercurrent is fluctuating arbitrarily and the state 
is characterized by a persistence time r for the super- 
currents: (J(O)J(t)) = ( J 2 ) exp — t/r. One may estimate 
that t is self-generated by the vortices themselves and is 
of order h/r = h 2 nv/m. The conductivity of such a vor- 
tex tangle will be a — p$T '. We may speculate that when 
h/r drops below kT, or equivalently when the number 
of vortices drops below a critical value where their en- 
tropy no longer compensates for their kinetic energy, the 
vortices evaporate: this is T c , described in a Kosterlitz- 
Thouless way as suggested by Leoi&. This provides a 
basis for the empirical rule proposed by Homesii, as well 
as for the observations of Timusk on anomalous increase 
of r in the pseudogap region-^. An even more specula- 
tive argument based on the vortex tangle can explain the 
Nernst observations. (Fig. [31 shows a heuristic first at- 
tempt at a description of the Nernst observations). The 
fact that T c is controlled by the vortex liquid transition 
invalidates most intuitions about it from BCS theory-for 
instance, it makes the d-wave T c insensitive to scattering. 
It is from the locking principle that the two insights men- 
tioned in the Introduction arise. Why the Cu02 planes? 
Because they have the feature that nearest-neighbor ex- 
change with only four neighbors allows the two almost 
degenerate gap functions of even and odd symmetry in 
x and y, of which one may be used to enhance the ki- 
netic energy, the remaining one giving a strong x 2 — y 2 
pairing energy. Of course, there are other aspects, par- 
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ticularly the Jahn- Teller distortion which enhances the 
energy scale, and the fact that Cu +++ does not self-trap, 
and all mean that unfortunately the scenario is unlikely 
to be repeated. Why the nodes? Because the RVB can 
only be a dyad: the spin interaction does not have a third 
possibility for pairing. Thus only one function can be left 
over as a gap function, and it must have nodal lines which 
do not lie along the Fermi surface. 



B. Hole-Particle Asymmetry 

One of the more significant experimental anomalies 
of the cuprates is the marked hole-particle asymmetry 
of the vacuum tunneling spectra. To those of us who 
worked on BCS superconductivity theory, this is par- 
ticularly striking because it is never observed in those 
materials. There is a large "peak-dip-hump" structure 
observed on the side on which holes are injected, be- 
coming stronger as the sample is underdoped (see Fig. 
0] ). The underlying band structure is not responsible 
since it is theoretically irrelevant and experimentally im- 
plausible. In tunneling, a theorem of Schrieffer removes 
much of the effects of quasiparticle interactions, so that 
the broad spectra seen in ARPES are referred back to 
the quasiparticle pole energies; the "hump" structure in 
fact has a strong resemblance to the incoherent part of 
the ARPES EDC's. It is a remarkable achievement of 
"plain vanilla" that it can give a sometimes quantitative 
account of these spectra. In order to do so we must mod- 
ify the ansatz Eq. 101 for <f> in Eq. (@J. BCS functions 
are wave-packets in the space of total electron number 
and one makes up non-number-conserving quasiparticles 
by taking advantage of this fact. This grand canonical 
approach is justified because the packet is centered at the 
correct particle number and the amplitudes for N — 2, N 
and N + 2 are essentially identical. But the projection 
process, while it does not change particle number, does 
project out very different numbers of states, so that after 
projection the wave packet is skewed in A-space. In or- 
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FIG. 3: A model calculation for the Nernst observations. 
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FIG. 4: Tunneling spectrum in optimally doped BSCCO. 
Data from S. H. Pan (unpublished). 



der to move the center of the packet back to N, we must 
introduce a fugacity factor dependent upon N: 
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-(«T+«i)/2 



(7) 



and g turns out to be the familiar kinetic energy renor- 
malization factor 2x/(l + x). Although in Q it is clear 
that the factor g cannot change any energy calculation 
since the Hamiltonian and projection conserve particle 
number, it is vital in understanding the process of tun- 
neling where a particle is added or removed. Eq. Q 
may be rewritten by distributing the factors of g among 
the terms of the product, appearing very different but 
actually this is an obvious identity: 



$ = Y[(u k + v k c^cl kl )\0) 

k 

<x Y[(u k +Wfc4 T cL fei )|0) 

k 



(8) 



where, u k = guk/ y/g 2 ul + v'f, and v k = v k /y/g 2 u 2 k + u|. 
In JHJ , the ratio of probabilities of zero and single occu- 
pancies is correct for the projected state and is thus not 
altered by projection. What it makes clear is that the 
projected state is made up from singlet pairs in which 
the relative amplitude of paired holes (the u term) is de- 
creased relative to that of paired electron spins (the v 
term) by the factor g. In a sense, there are two types 
of condensed bosons, the valence bonds of the RVB and 
the hole pairs, and in this theory we set their relative 
amplitudes free, although they remain coherent: they 
are "locked" together. The principle on which we cal- 
culate the tunneling spectrum is the following, Once we 
have chosen the form ||SJ| for we may define the single- 
particle excitations whose energies satisfy the gap equa- 
tions in terms of the wave functions, Pci CT $ and Pc icr $>, 
or equivalently, Pc k( j& and Pc k <f>, and these are now 
equivalently normalized. But the matrix elements of the 
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tunneling process insert a particle or a hole prior to the 
projection operation, at a particular site effectively, so 
that they connect to the operators Ci a P and c ia P, and 
we have to commute the fermion operator through the 
projection operator to determine its effect. 

We may write = Pc* + (1 — P)c) where (1 — P) 
projects onto states with a doubly-occupied site which 
are effectively at infinite energy (after the canonical 
transformation) . Thus whenever the inserted particle en- 
counters an occupied site, the state is projected out, and 
only with probability x does it encounter an empty site, 
i.e., Pc* . But when it is Pc^ it lands in a legitimate 
excitation, i.e., Pc^P — Pet. Thus when a particle en- 
ters (with probability x) it does so coherently. The hole 
problem is less obvious, c may be commuted through P 
with the result, 

= P ^ ~ AT 51 nk -° ~ H C l-a C k'~a)c la ■ (9) 
k k^k' 

The second term, when acting on is simply a number 
times c: (rii- a ) = (1 — x)/2. The third term is genuinely 
incoherent, creating three excitations; but these three can 
come from any energy in the spectrum so we expect this 
term to be quite small everywhere and to rise only as the 
square of the tunneling voltage for small voltages. The 
net effect of ijjjj, then, is that 

CP W ^~Y~PC ■ (10) 

Thus the ratio of the probabilities of tunneling of elec- 
trons vs holes is (no surprise!) g = 2x/(l + x). At high 
energies w> A, where the quasiparticles are pure holes 
or electrons, this is the expected asymmetry, and insofar 
as experiment is able to ascertain, apparently this ratio 
agrees well (taking into account the small error caused 
by the canonical transformation expiS). 

The spectrum at lower energies is complicated by the 
fact that superconducting quasiparticles are mixtures of 
electrons and holes. At the Fermi surface, exactly at the 
gap energy, they are equal mixtures and the singularity 
at the gap must be identical for the two sides. The work- 
ing out of the exact interpolation formula for the tunnel 
current is a little complicated and I give here only the 
formulas: the tunneling density of states for electrons is 

Ne{E , A) = ^-g( U " + V " ) . (11) 
dE y^u 2 + v 2 g 2 \/v 2 + u 2 g 2 J 

The g factor in this formula comes from the projection 
factor-which I emphasize does not multiply the matrix 
element, it is essentially a relative number of open chan- 
nels. As we see, for u k 1, at high voltage, the limiting 
value, 1, comes from the second factor and the tunnel- 
ing is suppressed by g. On the other hand, for holes the 



tunneling density is 

N h (E, A) = ±g ( f + f ) . (12) 

Here the g factor comes from the normalized fugacity 
factor, and at high voltage u « 1, v w and g cancels 
out, giving the ratio g between the two limits. These 
formulas fit data surprisingly well. 

These are the formulas for fixed A. Note that for e = 0, 
at the Fermi level, u = v and the two are identical, the 
"coherence factor" amounting to gj \J\ + g 2 . This agrees 
with sum rule arguments. The asymmetry begins, how- 
ever, with a vertical slope at A, so cannot to be said to 
be exclusively a background phenomenon, as is seen most 
clearly in the fact that the peaks of observed spectra (see 
Fig. appear to sit on background levels of different 
heights. These formulas must be integrated over the d- 
wave distribution of gap values to give a prediction for 
comparison with observed spectra. This we have done 
only roughly, using P(A) = — A 2 , as though the 

Fermi surface were circular and not taking into account 
the actual band structure, which does somewhat affect 
the distribution of A values. In Fig. [B] we give the pre- 
dicted spectra for a number of values of <?, using this 
simplification. 

The fit to experiment, at least in the main features, 
is fundamentally significant. Of course, it helps con- 
firm the basic structure of the theory, and the use of 
superexchange as the major pairing interaction. But it 
has even deeper implications. One is that even though 
it is basically a mean field theory based on an Hartree- 
Fock ansatz, it is not a Fermi liquid-based theory, that 
is to say that in no way can it be adiabatically contin- 
ued to a BCS-like modification of Fermi liquid theory. 
The most fundamental property of Fermi liquid theory 
is hole-particle symmetry«2i; after all, how can one have 
a theory based on a distribution of quasiparticles unless 
that distribution counts particles minus holes, 1 for 1? 
This projective theory has destroyed that symmetry in a 
very fundamental-yet simple- way. 

Yet the projective feature is rooted in the real physics 
of the system. As pointed out also by Capello et at, 
once one is above the Mott critical U c , there are what 
we called "anti-bound states'— and what Ref. ^01 calls 
"holon-doublon bound states" which cannot be treated 
perturbatively but must simply be projected out of the 
problem. A result is that the spectrum is overcomplete 
and the particle operators do not obey a simple fermion 
algebra. I believe that the broad features seen in the 
momentum-resolved particle spectra are related to that 
problem, but fortunately the tunneling spectrum is sim- 
plified by the benign results of "Schrieffer's theorem" and 
is easier to interpret. 
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ergy of the pairs; or it is possible that there is a tendency 
to bistability near the low-doping quantum critical point. 

Finally, there is the question of phase fluctuations of g 
(which is the appropriate variable to assign a phase to, 
since it controls the charge carriers.) It is known that 
the phase transition at T c is of "A — Y" character both 
for optimal dopingi^ and for very low doping near the 
quantum critical points^ , and that above T c there is 
a large region in which the state is best described as a 
"vortex liquid" rather than a normal metal, i.e., there is 
a fluctuating superconducting order parameter, (see Ref. 
fl5h This has been described as a regime in which £ is still 
locked to the kinetic energy but A, i.e., the phase of g, 
is freely fluctuatingii. There is a very important open 
question here as to whether or not there is a transition 
into a still higher T phase which has an RVB but is not 
a vortex liquid^. 



IV. DISCUSSION: ALTERNATIVE 
APPROACHES 



FIG. 5: Predicted tunneling spectra for various Z(= g). 



C. Fluctuations of the Asymmetry Parameter g 

The asymmetry parameter, as I remarked, plays a role 
similar to that of a condensate of hole pair bosons, which 
is locked to the spinon pairs of the RVB by the charge- 
spin locking process. Its average value is determined by 
the doping and charge neutrality, but it is evident that 
formally we can allow it to vary either in space or in 
momentum space. More speculatively, we can allow it a 
dynamic character, and I believe that its collective modes 
do in fact give us extra degrees of freedom which play a 
role in the non-Fermi liquid behavior mentioned above. 
This and the other remarks I will make in this rather 
speculative section are inspired more by suggestions from 
experimental observations than by apriori theory, but I 
do believe that treating g as a physical object can lead 
to considerable insights. 

The most obvious is the possibility of allowing g to 
vary along the Fermi surface, just by making it a func- 
tion of k in (JSJ). The mean value x = (<?/(2 — g)) must 
be maintained for charge neutrality. The phenomenon 
of "Fermi arcs" is observed by ARPES in underdoped 
systems, where the regions of the Fermi surface near the 
nodes remain sharply defined while the antinodal regions 
smear out and disappear. The nodal regions are also 
those where the kinetic energy is greatest, so that one 
could gain energy by making g large at the nodes and 
small at the zone corners. No calculation of this pro- 
posed effect yet exists. 

It has also been proposed that the hole percentage may 
vary spatially, in particular that at low doping g could 
form a kind of charge density wave or "superconducting 
electron solid"—. The motivation could be Madelung en- 



The RMFT theory works, based on t — J physics and 
superexchange as the interaction, and can account semi- 
quantitatively for the basic phenomena of cuprate su- 
perconductivity, and qualitatively for many more. Why 
then are contradictory theories being promoted? The 
most popular theories reject the Mott- Anderson physics 
entirely and go in contrary directions. There seems to 
be a psycho-social need among physicists for an explana- 
tory boson, some kind of tangible glue to hold the pairs 
together, I suppose because of the folk memory of talks 
about the BCS mechanism and the analogy of two bod- 
ies on a mattress; or else a simplified view of Feynman 
diagrams. It is felt, I suppose, that the Mott theory 
is based on purely repulsive forces-but those of us who 
actually worked on BCS recognize that the phonon inter- 
action is not literally an attraction either, merely a par- 
tial screening of the electrons' Coulomb repulsion. Why 
a superexchange integral universally agreed and experi- 
mentally measured to be of order 1000 degrees is thought 
to be inadequate for pairing has always escaped me; but 
it is. The two most popular glues are phonons and anti- 
ferromagnetic spin fluctuations. 

Phonons start out with a big disadvantage: the BCS 
concept is irrevocably based on an on-site, local interac- 
tion; and is incompatible with d-wave. In the cuprates, 
the phonons are undoubtedly optical ones involving the 
oxygen octahedron (oh, there are other suggestions, but 
even less plausible) and there are perhaps ways of distort- 
ing these in order to give a d-wave, but I have never seen 
a plausible one. Intrinsically, Einstein optical phonons 
lead to local interactions. But, experiment is the best 
teacher. The isotope effect measurements of Keller— find 
a reasonably-sized isotope effect on T c , apparently con- 
firming the phonon hypothesis; but Keller was thorough 
enough to also measure the isotope shift of p s , the su- 
perfluid density, cx A~ 2 ; and he finds that this shifts by 
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the same fractional amount. It was pointed out very 
early in the game by Fisher et al.jH that unlike the poly- 
electronic metals for which BCS theory works and the 
isotope shift comes entirely from the pairing interaction, 
oxides are best understood as tight-binding systems with 
interactions which depend exponentially on interatomic 
distances. Thus zero-point vibrations will have an appre- 
ciable effect on normal state properties such as the band 
mass which determines p s . Since T c is an X — Y tran- 
sition as already remarked, its value is expected to be 
directly proportional to p s , the coupling in the xX — Y 
model, as observed, so that apparently there is no experi- 
mental isotope shift ascribable to the pairing interaction. 
In fact, even if there were, T c is insensitive to the actual 
value of A, as we explained above. Extensive ARPES 
studies have catalogued what may be phonon effects on 
the quasiparticle dispersion^ but these seem to be irrele- 
vant to the pairing mechanism. J, of course, itself varies 
in a similar way as t with interatomic distance and may 
provide a partial source of the observed isotope shifts 
in dispersions. It seems that calculating phonon effects, 
while worth doing for its own sake, is not the most urgent 
task. 

There are other phonon schemes, most notoriously 
the bipolaron theory. One understands the impulse to 
look this way, since polaron phenomena are so ubiqui- 
tous in oxides. But very early on it became clear that 
one reason the cuprates are so favored is that this case 
is gloriously free of polaron effects, presumably because 
Cu ++ and Cu +++ have similar Jahn- Teller displace- 
ments. The remarkably detailed tunneling and ARPES 
spectra demonstrating well-characterized quasiparticles 
exclude small polaron phenomenology. I believe that 
Baskaran's theorji2£ explaining the electron-doped case 
as dominated by small polarons must be essentially cor- 
rect, and the contrast with hole doping illustrates well 
what phenomenology polarons might lead to. 

A second putative source of the "glue" boson is "anti- 
ferromagnetic spin fluctuations" . This idea sounds simi- 
lar to the Mott-based theory but is not at all so, in fact 
proceeds on exactly the opposite principle: that in the 
end the physics is to be obtained by "summing all the 
diagrams" starting from a Fermi liquid^. Another way 
to say it is that the assumption is that the theory fits 
under the general scheme of Ref. [§, where all interac- 
tion terms are renormalized downwards, while the plain 
vanilla theory makes the assumption that one must start 
by renormalizing U — ► oo, with the Rice canonical trans- 
formation. I feel that U c , the Mott critical U, marks 
a fundamental separatrix between basins of attraction, 
and that the cuprate case is on the large U side. The key 
question is whether the frequency associated with most of 
the pairing interaction is above a Mott-Hubbard gap, and 
therefore cannot be represented by a boson whose spec- 
trum extends continuously to zero frequency. In that case 
it might as well be represented by a simple four-Fermion 
vertex J. The idea of antiferromagnetic spin fluctuations 
is that the opposite is the case, and that somehow if 



one can sum enough diagrams the Mott gap will disap- 
pear from the problem and interactions will proceed by 
the exchange of a putative low-frequency spin-fluctuation 
boson. 

Since, in fact, one cannot come close to summing all 
the diagrams, papers based on this idea have tended 
to contain about one parameter per experimental fact, 
and therefore to "explain" great numbers of these facts. 
Apparently recent advances in experimental detail have 
led to exhaustion of invention, and many rather cru- 
cial discoveries remain unexplained, for example Fermi 
arcs, tunneling asymmetry, the vortex liquid phase, the 
checkerboard, Homes' identity. 

There are a number of more mysterious suggested 
sources for the "glue boson" , many of which invoke the 
equally opaque concept of a hidden "quantum critical 
point" ; their variety excludes detailed explication. 

Perhaps no longer worthy of mention is the "stripe 
theory" , the problem of which was that it never seemed 
to be a theory of the superconductivity, but only a theory 
of the stripes themselves. Since stripes are not common 
to many of the cuprate superconductors, and as time 
goes on to fewer and fewer, it is hard to understand their 
relevance. 



V. CONCLUSIONS AND ANTICIPATIONS 

Many of my conclusions were rather strongly stated in 
the Introduction. It seems that the Gutzwiller method 
works perhaps even better than we had any right to ex- 
pect. It also has the added feature that it brings out 
the deep difference in principle between a Fermi-liquid 
based approach and the actual behavior of the cuprates 
in a relatively simple and straightforward way, both in 
demonstating the hole-particle asymmetry of the Green's 
functions and in the "locking" phenomenon. 

Quite understandably, there are other ways to ap- 
proach the same physical model, and some of them have 
a good chance of being more accurate or rigorous-for in- 
stance gauge and slave boson theories, one of which I 
quoted here. One can certainly differ on the applicabil- 
ity of the crude approximations made in Plain Vanilla 
to make it soluble; and it is very meaningful to try to 
add in further terms to the interactions used, and to 
study the various accompanying phenomena such as co- 
existence with antifcrromagnetism. The major puzzle re- 
mains that of the Strange Metal, the mysterious phase 
above T* , and the strange quantum critical point where 
the d-wave gap goes to zero. The linear T, linear in 
electron-electron scattering mechanism which pervades 
the high-energy region is still a puzzle but must be a 
characteristic of the purest Mott physics 32 . 

I should not fail to mention the accumulation of recent 
direct or semi-direct calculational results all of which are 
now tending to converge on the conclusion that d-wave 
superconductivity undoubtedly appears in the Hubbard 
and t — J models. I am sure these will be represented well 
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elsewhere in this volume; and of course, I have absolutely 
no problem with them; it is a matter of taste whether one 
prefers approximations such as Plain Vanilla which allow 
understanding of the phenomenology, or more exact but 
only semitransparent calculations. 

Throughout the paper I have alluded to avenues for 
further exploitation of the method, specifically the possi- 
ble explanation of the "Fermi arcs" as a /c-dependence of 
g, and of nanoscale structures as spacial modulations of 
it; but both will require more detailed calculations than 
we are yet capable of. 
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